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ABSTRACT 

The  aim  of  this  project  is  to  develop  an  understanding 
of  the  mathematics  and  physics  involved  in  a  bungee 
jump,  in  order  to  gain  a  better  understanding  of  the 
concepts  involved.  We  will  use  Newton’s  law  and 
utilize  techniques  needed  to  solve  second  order 
differential  equations.  For  simplicity;  we  will  regard 
the  jumper  as  a  point  mass  in  one  dimension  and  we 
will  make  the  assumption  that  wind  resistance  has  a 
linear  relationship  to  speed.  We  will  also  assume  that 
the  bungee  cord  obeys  Hooke’s  laws  and  doesn’t 
break. 

Keywords:  Newton’s  second  law;  Hooke’s  law; 
Young’s  modulus;  Wind  resistance;  Driving  force; 
Terminal  velocity. 

Modeling  the  Motion 

With  a  bungee  jump,  the  model  of  the  motion  in  one 
dimension  with  respect  to  time  can  be  broken  into  two 
different  sections:  the  free  fall  from  the  ledge  and  the 
time  when  the  bungee  cord  is  pulling  on  the  jumper. 
To  model  both  sections  we  need  Newton's  second  law 
of  motion.  It  says  that  the  total  of  all  the  forces  acting 
on  an  object  is  equal  the  mass  of  the  object  multiplied 
by  its  acceleration. 

TF  =  ma 

To  describe  how  the  bungee  cord  will  pull  on  a 
bungee  jumper  we  will  use  Hooke's  law.Hooke's  law 
says  that  the  force  exerted  by  a  spring  is  proportional 
to  the  distance  that  it  has  been  stretched  from  its 
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equilibrium  point.  The  relationship  is  said  to  be  F  =  - 
kx,  where  k  is  a  constant  that  is  specific  to  the  spring. 
The  negative  is  in  the  equation  so  that  the  force  is  in 
opposite  direction  of  the  stretch  (pushing/pulling 
toward  the  resting  point  of  the  spring).  For  our  model, 
we  assume  that  a  bungee  cord  will  behave  like  aspring 
described  by  Hooke's  law.  The  spring  constant  k  can 
be  found  experimentally  by  hanging  a  known  weight 
on  a  spring  and  measuring  how  far  it  stretches,  also,  k 
can  be  found  by  rearranging  the  Young's  modulus 
formula.  The  Young's  modulus  formula  applies  to 
objects  that  have  a  constant  cross-sectional  area  (same 
thickness  throughout)  and  objects  that  have  the  same 
material  all  throughout,  the  formula  also  only  applies 
for  forces  that  compress  or  stretch  the  object.  The 
formula  says  that  any  material  has  constant  that 
relates  the  stress  on  an  object  to  the  strain  on  the 
object.  The  stress  is  said  to  be  the  force  applied  on  the 
object  divided  by  the  cross-sectional  area  of  the 
object,  F/A.  The  strain  is  said  to  be  the  change  in 
length  divided  by  the  initial  length,  dL/L.  The 
Young's  modulus  formula  is: 

y-f/a 

dL/L 

where  Y  is  the  Young's  modulus,  a  constant  for  a 
given  material.  F  is  the  force  being  applied,  A  is  the 
cross-sectional  area.  L  is  the  length  of  unstretched 
cord,  and  dL  is  the  change  in  the  length  of  the  cord 
once  it  has  been  stretched.  This  can  be  rearranged  to 
look  like  Hooke's  law,  if  we  remember  that  dL  and  the 
x  from  hook's  law  have  the  same  meaning. 
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The  Young's  modulus  constant,  the  cross-sectional 
area,  and  unstretched  cord  length  will  not  change  for 
any  given  cord,  so  we  can  say  that  the  spring  constant 
in  Hooke's  law  is: 


k 


d2x  B  dx 

-  „  H - —  =  0 

dtz  m  dt 

The  characteristic  equation  for  the  homogeneous 
equation  is: 


p  (A)  =  A2  +  —  A 

m 


In  order  to  keep  the  math  doable,  we  will  also  stick 
with  a  wind  resistance  that  is  linearly  proportional  to 
the  velocity.  There  are  more  accurate  and  complicated 
formulas  to  account  for  wind  resistance,  but  for  this 
paper  we  will  be  using: 

Fwind  —  —fiv 

In  this  paper,  we  will  define  the  positive  x  direction  to 
be  going  down. 

Free  fall 

When  the  bungee  jumper  first  jumps  off  of  the  bridge, 
they  will  have  the  forces  of  gravity  and  wind 
resistance  acting  on  them.  The  force  of  gravity  is 
positive  because  we've  defined  the  positive  x  direction 
to  be  going  down.  The  wind  resistance  is  —  /?vbecause 
the  force  of  wind  resistance  is  always  in  the  opposite 
direction  of  the  velocity  and  because  we  decided  that 
making  the  force  of  wind  resistance  proportional  to 
velocity  was  appropriate  for  this  paper.  Using 
Newton's  second  law: 

mg— pv  =  ma 

using  the  fact  that  a  =  dv/dt  ;  v  =  dx/dt  and  a  little 
algebra,  we  can  put  this  into  the  form  of  an 
inhomogeneous  equation. 


d2x  P  dx 
dt2  m  dt  ® 


(2) 


Solving  the  homogeneous  equation 

To  find  the  general  solution  of  the  inhomogeneous 
equation,  we  must  first  find  the  solutions  to  the 
equivalent  homogeneous  equation.  In  other  words,  we 
will  set  the  driving  force  g  to  zero.  (Note:  there  are 
other  easier  ways,  but  we  will  solve  our 
inhomogeneous  equation  this  way  so  that  the  process 
is  the  same  as  the  process  needed  in  the  next  section.) 


We  set  p(A)  to  0  to  find  that  the  roots  of  the  equation 
are: 


The  solution  to  homogeneous  equation  is: 

xh  =  C1eAlt  +  C2eXlt 
Substituting  in  our  values  for 

_££ 

Xh  =  C1  +  C2e  rn 

Solving  for  the  particular  solution 

Now  that  the  homogeneous  equation  has  been  solved, 
the  only  thing  left  to  find  is  the  particular  solution.  If 
we  look  at  equation  (2)  and  we  look  at  the  driving 
force  that  is  on  the  right-hand  side,  then  we  see  that 
the  driving  force  is  a  polynomial  of  zero'thdegree  (a 
constant).  The  particular  solution  to  the  equation 
should  be  a  polynomial.  If  we  say  a  particular  solution 
is: 

xv  =  at  +  b  ,—  xv  =  a  and  ^xv  =  0 

p  ’ at  p  dt2  p 

Then  we  can  plug  this  into  (2)  and  solve  for  a 

(0)  +  — (a)  =  g 
m 

p , , 

—  (a)  =  g 


Our  general  solution  is: 

r  i  n  i  m3  *.  ,  u 

X  =  C,  +  C7e  rn  +  —  t+b 

J.  z  p 
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Ct andb  are  both  constants,  because Cj  is  an  arbitrary 
constant,  we  say  that  it 
accounts  for  b. 


pt 

x  =  +  C2e  m  + 


mg_ 

P 


(3) 


{P/m)2  —  4 (/c/m) is  greater  than  zero,  equal  to  zero, 
or  less  than  zero. 

Case  1 

{p/m)2  -  4  (^)  >  0 


Below  the  equilibrium  point  for  the  cord 

So  far,  an  equation  for  the  motion  of  the  bungee 
jumper  has  been  found.  Now,  we  need  to  find  an 
equation  of  motion  for  when  the  bungee  jumper  is 
below  the  equilibrium  point.  To  do  this  we  look  back 
to  Newton’s  second  law: 

£F=ma 


This  is  said  to  be  over  damped,  and  the  solution  will 
not  oscillate.  If  bungee  cord 

and  jumper  system  fall  into  this  category  then  there  is 
a  problem,  because  the  bungee  jumper  will  likely 
experience  too  much  acceleration.  In  this  case,  the 
characteristic  equation  has  two  distinct  real  roots.  The 
solution  will  take  the  form, 


When  the  jumper  is  below  the  equilibrium  point,  they 
still  have  gravity  and  wind  resistance  acting  on  them. 
But,  now  they  also  have  the  force  of  the  cord  pulling 
on  them.  We  will  say  that  the  force  applied  by  the 
bungee  cordis-kx  (from  Hooke’s  law).  Using  this 
with  Newton’s  second  law: 


„  dx  ,  d2x 

mg-p--kx  =  m  — 


Once  again,  this  need  to  be  rearranged  algebraically  to 
look  like  an  inhomogeneous  equation 

d2x  p  dx  k 

tt- I - 7-  H — x  =  q 

dt *  mdt  m 


Solving  the  homogeneous  equation 

The  first  step  to  solving  the  inhomogeneous  equation 
is  to  solve  the  equivalent  homogeneous  equation. 


d2x 

dt2 


P  dx  k 

-\ - —  -I — x  =  0 

m  dt  m 


The  characteristic  equation  is 


xh  =  C \ex^  +  C2eA2f 
We  can  plug  in  our  expression  for  A  to  get, 

(-/?/jn±J(/?/m)2-4(fc/m))t/2 

xh=C1e  't _  + 

(-p/m±l{p/m)2-4(k/m))t/2 

C2e  ’ 

Case  2 

{P/m)2  —  4  {k/m)  =  0 

This  case  is  said  to  be  critically  damped,  and  its 
solution  appears  very  similar  to  the  over  damped 
case.The  characteristic  equation  has  one  real  root,  and 
the  solution  will  have  the  form. 

xh  =  C1eAt  +  C2teXt 
Or 

Xn  =  (Ci  +  C2t)eAt 
Remembering  that  in  this  case 

{P/m)2  —  4 {k/m)  =  0,  we  say  that  A  =  —p/2m  and 
thatour  solution  is, 


P(A)  —  A2  +  —  A  H — 

mm 


xh  =  {C,  +  C2t)e~^2m 


We  set  P(A)  to  0  and  use  the  quadratic  to  find  that  the 
roots  of  the  characteristic  equation  are: 

^_-P/m±i/  QS/?n)2-4(fc/m) 


At  this  point,  the  solution  to  the  homogeneous 
equation  can  go  in  three  directions,  depending  on  the 
discriminant.  There  are  three  cases,  when 


Case  3 

{p/m)2  -  4  (^)  <  0 

This  case  is  said  to  be  under  damped,  and  it's  solution 
oscillates.  The  solution  in  this  case  will  be  a 
combination  of  the  real  and  imaginary  parts  of  eAt.  To 
find  this,  we  need  to  know  what  the  real  and 
imaginary  parts  of  e(a+lb)t  are.  To  do  this  we  first 
split  up  the  exponent. 
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e(a+i  )t  _  eate(bi)t 

efcltsplits  up  into  cos  bt  +  i  sin  bt,  so  we  get: 

eate(bi)t  _  gat^cos  bt  +  i  sin  bt) 

=  eat  cos  bt  +  ieat  sin  bt) 

The  real  part  of  the  equation  above  is  eatcosbt  while 
the  imaginary  part  is  eatsin  bt.  If  plug  in  the  A  value 
that  we  got  from  equation  (4),  and  we  use  the  same 
process  to  separate  out  the  real  and  imaginary  part, 
then  we  get  that  the  real  part  of  eh  is: 


Because  there  is  no  t  term  the  coefficient  a  must  be 
zero  and  as  a  result: 

k 

—  b  =  g 
m 

b=^ 

k 

So  the  particular  solution  is: 

gm 

x”=— 


e-nt/2mcostMMEhiiM 

While  the  imaginary  part  is: 

_Btnm  .  J\(J3/m)*-4(k/m)\ 

e  pc/zmsin  t - - - 

2 

The  homogeneous  solution  will  be  a  combination  of 
the  real  and  imaginary  parts  of  eXt 
xh 

_JL 

=  Cxe  2 m  COS  t 

+  C2e  2 m  sin  t 

Solving  for  the  particular  solution 

The  next  step  is  to  get  the  particular  solution  of  the 
inhomogeneous  equation. 

d2x  (3  dx  k 

-rr~ I - T  H — x  =  g 

dtz  mdt  m 


PROBLEMS 

1.  Solving  without  a  numerical  solverSetting 
variables 

Before  we  can  take  the  last  steps  to  solving  for  the 
position  as  a  function  of  time,  we  need  some  initial 
conditions  and  parameters.  We  are  going  to  model  an 
80kg  man  jumping  of  a  bridge  with  a  height  of  100m 
with  a  30  meter  long  bungee  cord  with  adiameter  of 
.00635  meters  (quarter  inch). 

We  measured  a  bungee  cord  that  had  the  same 
thickness  and  a  length  of  1.31m  and  found  that  when 
a  force  of  10N  was  applied  the  cord  stretched  .16m. 
The  bungee  cordcould  stretch  to  a  maximum  of  220 

SOLUTION 

Y  -  F/A 
dL/L 

Y  _  10(tt(0.00635/2)2 
0.16/1.31 

=  2.59  x  106 


Vl(£M)2-4(/c/m)| 
_ 2 _ 

j\(J3/m)2  -  4(k/m)\ 

2 


Fortunately,  the  particular  solution  does  not  have 
multiple  cases  to  consider.  In  addition,  the  driving 
force  is  a  polynomial  of  zero'th  degree  just  like  it  was 
for  the  free  fall  equation.  We  can  assume  that  the 
solution  is  some  polynomial  with  unknown 
coefficients  and  then  plug  the  solution  into  the 
equation  above  to  find  those  coefficients. 

d  d2 


xp  =  ax  +  b, 


=  a, 


=  0 


We  will  assume  that  the  cord  that  we  measured  is 
suitable  for  bungee  jumping.  We  will  also  assume  that 
they  could  will  stretch  linearly  with  the  force  applied 
so  that  we  can  use  the  Young's  modulus  as  well  as 
Hooke's  law.  With  these  assumptions  the  k  value  that 
is  used  in  Hooke's  law  will  now  be  determined. 
Pulling  from  equation  (1): 


K  = 


AY 

L 


plugging  these  into  the  inhomogeneous  equation  we 
get, 

(0)  +  —  a  +  —  (ax  +  b)  =  g 

m  m 


_  (7t(0.00635)2)2.59x106 

Jv - 

30 

K  =  2.73 


@  IJTSRD  |  AvailableOnline@www.ijtsrd.coml  Volume -2  |  Issue -2  |Jan-Feb2018 


Page:  867 


International  Journal  of  Trend  in  Scientific  Research  and  Development  (IJTSRD)  ISSN:  2456-6470 
Now  we  have  all  but  one  variable,  /?.  To  get  this  Let's  replace  -p/m  with  A  to  make  the  next  few  steps 


variable  we  can  solve  for  based  on  accepted  values  for 
the  terminal  velocity  of  a  person.  To  do  this,  we  go 
back  to 

2  F  =  ma  and  look  at  free  fall. 

easier 

x(t)  =  Ci  +  C2  e  ^  H — -jj- 1 

n  dv 

mg-pv  =  m- 

Differentiating  we  get, 

Because  we  looking  at  when  the  person  has  reached 
terminal  velocity  (no  change  invelocity),  we  can  set 
dv=dt  to  zero  and  solve  for  /?. 

d  „  mq 

—  x(t)  =  v(t)  =  —C2Ae  +-jf 

If  we  use  the  initial  conditions 

mg  -  Pv  =  0 

x(0)  =C1  +  C2e^(0)+^(0) 

p  =  !M 

V 

v(0  )  =  -C2Ae~AW  +^ 

In  solving  for  p  we  had  already  assumed  that  the 
person  was  at  terminal  velocity,  so  we  can  plug  in  the 
value  that  we  looked  up  for  v.  According  to 
hypertextbook.com,  a  good  approximate  of  the 
terminal  velocity  of  a  person  is  60m/s.  with  this  we 
get: 

This  simplifies  to, 

x  (0)  =  Cl  +  C2 
v(0  )J-AC2  +  rf 

_  (80)(9.81) 

P  60 

Solving  for  C1  and  C2 

p  =  13.08 

Ci  =  x(0)  —  C2 

RESULT 

_  v(0)-(mgy<J3)) 

2  A 

YOUNG’S  MODULUS  (Y)  =  2.59  x  106 

HOOK’S  CONSTANT  (K)  =  2.73 

VARIABLE  (p)  =13.08 

m(v(0)  -  (mg)/ (/?)) 

Ci  =  x(0) + 

2.  At  the  start  of  the  jump 

n  m(v( 0)  -  (mg)/(p)) 

C2~  p 

Now  that  we  have  all  of  the  variables,  we  can  use  the 
use  the  general  solutions  that  we  found  earlier.  First 
we  will  use  the  initial  conditions  of  x(0)  =  -30  (30m 
above  the  equilibrium  point  of  the  unstretched  cord) 
and  v(0)  =  0.  Since  we  are  using  x(0)  and  y(0),  the 
initial  t  value  has  to  be  zero.  We  have  the  equation  (3) 
with  the  unknowns  C1  and  C2.  To  solve  for  both 
unknowns  we  will  have  to  use  equation  (3)  and  its 
derivative. 

If  we  plug  in  the  known  values  for  our  parameters  and 
conditions 

_  ,  80(0  -  (80)(9.81)/ (13.08)) 

Cl"  30  +  13.08 

r  80(0-(80)(9.81)/(13.08)) 

Ly - 

z  13.08 

SOLUTION 

From  this  we  get  that  Cx  =  -396:97  and  that  C2  = 
366:97.  So,  the  equation  for  free  fall  becomes 

x(t)  =  C1  +  C2e  m+-^t 

x(t)  =  366.97 e-°  1635t  +  60 1  -  396.97 
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From  this  position  function  we  need  two  things.  First, 
we  need  to  find  when  x  is  greater  than  zero  because 
this  is  when  we  need  a  new  equation  to  represent  the 
bungee  jumper.  Second,  we  need  to  find  the 
acceleration  as  a  function  of  time.  From  the  graph  of 
the  position  function  it  seems  that  x(t)  =  0  when  t  = 
2.67.  If  we  plug  this  and  all  of  our  knowns  into  the 
velocity  function: 


v(t) 


—C2Ae~M 


Then  we  get  that  v(2.67)  =  21.22.  We  will  need  this 
velocity  as  well  as  t  =  2.67  for  the  initial  conditions  of 
the  next  equation. 

RESULT 

TIME  (t)  =  2.67 
VELOCITY  v(2.67)  =  21.22 

3.  After  the  cord  starts  to  pull 

To  solve  this  part  of  the  situation  we  still  use  the  same 
parameters,  but  different  initial  conditions  and 
equations.  To  make  the  solution  easier  to  find  we  will 
say  that  the  time  2.67s  for  the  previous  equation  is  at 
Os  for  the  new  equation.  The  position  and  velocity 
stay  unchanged,  however.  For  this  part  of  problem, 
we  need  to  decide  which  case  to  use.  If  we  plug 
known  values  into  {P/m)2  —  4 (k/m)  then  we  get  - 
0.1098  which  means  that  we  must  use  case  3.  Our 


Now  we  can  put  this  into  a  matrix  and  row  reduce 

e~Ati  cos(5tj) 

(— Ae~Ati  cos(Rtj)  —  Be~Ati  sin  ^ Btt  e  A"  j) 


solution  for  x  is  a  combination  of  both  the  solution  in 
case  three  and  the  particular  solution. 

SOLUTION 


x(t)= 

1  2 _ 

r  ■  4-  VlO 3M)2-4(fcM)|  gm 

C2e  zm  sint  — - — — -  +  — 

*  2  k 


This  is  a  lot  easier  to  deal  with  if  we  let 
A  =  /? /2m  and 

B  =  (sqrt  |  (p /m)2  -  4(k/m)| )  /  2 
x(t )=C1e~At  cos (Bt)+C2e~At  sin (Bt)  +  pp 
We  differentiate  this  and  get  that: 
v(t)=C1(— Ae~At  cos(Bt)  — 

Be  —A 7s  i  n  Bt) + C 2  (—A  e  —A  tsmBt+  Be  —A  /cos  Bt) 

If  we  put  x  and  v  in  terms  of  an  initial  condition  then 
we  get 

x(ti)=C1e_i4ti  cos(Stj)  +  C2e~Ati  sin(Rtj) 

v(ti)=C1(— Ae~Mi  cos (Btt)  — 

Be  —A t is  i  n  Bti) + C2(— A e  —A  t is  i  n  Bti+ 
Be—AticosBti) 


e  Ati  sin  (Bt/)  x(te) 

(— Ae~Ati  sin  (Bt/)  +  Be~Ati  cos(Btj))  v  (tj) 


From  this  we  get  that  C1  =  -315:94  and  that  C2  =  - 
168.61.  The  position  function  after  the  bungee  jumper 
has  gone  past  the  equilibrium  point  of  the  cord  is 
x(t) 

315. 94e-016  cos(0.1636t)  — 

168. 61e-0  1635  sin(0.1636t)  +  287.5603 

Alternatively,  we  could  have  made  solving  for  the 
initial  conditions  on  the  second  equation  easier  by 
saying  it  started  at  time  zero.  The  constants  C1  and  C2 
would  turn  out  different,  and  the  t  would  have  to  be 
replaced  with  t  -2.67.  The  equation  would  still  have 


the  same  time  range.  This  turns  out  to  be  easier 
because  all  of  exponential  terms  go  to  1 ,  all  of  the  sine 
terms  go  to  zero,  and  all  of  the  cosine  terms  go  to  1. 

RESULT 

Cx=  -315.94 

C2  =  -168.61 
x(t)= 

315.94e_0  1635t  cos(0.1636t)  - 
168.61e-016  sin(0.1636t)  +  287.5603 
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4.  A  cord  s  20m  long  un-stretched  if  its  elasticity  is 
50%  and  the  spring  constant  is  20N/m,  what  is  the 
maximum  force  that  can  he  held? 

SOLUTION 

Use  the  formula  F  =  kx 

X  is  calculated  by  taking  50%  of  theun-stretched  cord. 
Therefore  x=10m. 

simply  plug  the  values  into  the  formula. 

F=  (20N/m)(10m) 

F=200N 

RESULT 

This  is  the  maximum  force  that  the  cord  can  hold  is 
200N. 

Using  a  Numerical  Solver  to  find  the  Motion  of  the 
Jumper 

In  order  to  get  the  correct  motion  of  the  jumper,  we 
need  to  find  a  way  to  model  the  position  as  they  move 
up  and  down  due  to  the  bungee  cord  after  they  leap. 
Using  the  origin  as  the  equilibrium  point  of  the 
weighted  bungee,  when  x  >  0  the  motion  of  the 
jumper  can  be  modeled  as  the  sum  of  the  forces 
equals  the  mass  multiplied  by  the  acceleration.  To 
model  the  oscillation,  we  will  treat  the  bungee  cord  as 
a  spring,  and  employ  Hooke's  Law  to  model  the  force 
of  the  bungee  on  the  jumper.  Hooke's  Law  is  used  for 
calculating  the  force  imparted  by  a  spring  on  an 
attached  object  and  is  defined  as  F=kx  .where  x  is  the 
displacement  from  the  force.  The  spring  constant  k  is 
given  by  dividing  the  force  in  Newtons  by  the 
displacement  caused  by  that  amount  of  force. 

There  are  three  forces  acting  on  the  jumper  when  they 
are  under  the  influence  of  the  bungee.  First,  is  the 
force  of  gravity,  which  is  simply  mg.  Second,  there  is 
the  force  due  to  the  bungee,  which  is  modeled  after  a 
spring  with  a  force  F  =  -  kx  with  k  being  the  spring 
constant  of  the  bungee  cord.  Third,  the  force  of  wind 
resistance,  with  a  wind  resistance  constant  /?,  whose 
value  is  determined  above  in  equation  (??),  which  has 
a  linear  relationship  with  the  velocity  of  the  jumper. 
The  calculations  involved  in  finding  the  true  value  for 
P  are  rather  complicated  and  beyond  the  scope  of  our 
understanding,  so  we  will  use  the  calculated  constant. 
Adding  up  the  forces,  we  get 


Where  -  kx  is  the  force  of  the  bungee  on  the  jumper, 
and  Pv  is  the  force  of  the  wind  resistance.  It  makes 
sense  that  as  the  magnitudes  of  kx  and  Pv  increase, 
they  will  eventually  overcome  the  force  of  gravity. 
When  this  happens  the  jumper  will  be  at  the  end  of 
the  downward  motion  of  the  oscillation,  and  will  be 
shot  back  up  from  the  force  of  the  bungee. 

Finding  the  Equation  of  Motion 

In  order  to  find  an  equation  of  the  motion  of  the 
jumper,  the  equation  above  must  be  broken  down  to 
components  that  make  up  each  force. 

mg-kx-Pv=  ma 

Moving  the  terms  around,  we  get  the  expression 
(3  k 

x"  -\ - x'  -I - x  =  g 

m  m 

This  gives  a  differential  equation  with  a  forcing  term 
g,  which  can  be  solved  with 

relative  ease.  In  order  to  utilize  MATLAB's  pplane8 
utility,  we  must  have  two  first  order  differential 
equations  to  relate  to  one  another.  These  equations 
were  found  using  the  following  steps: 
substituting  a  =  (3/m  and  b  =  k/m,  and  isolating  the 
x"  we  get 

x"  =  g  —  ax'  —  bx 

Now,  the  next  step  is  to  create  two  first  order 
equations. 

Vi  =  x 
y2  =  x' 

Now  we  plug  yx  and  y2  into  our  equation,  with  the 
addition  of  y[  =  y2 

y'i  =  T2 

y'2=  g- ay 2-  by± 

When  a  and  b  are  substituted  back  into  the  equation 
we  get, 

y[  =  Y2 

,  P  k 
y2  =  g  y2  y± 


r gravity 


—  kx  —  (3v  =  F1 


net 
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RESULT: 

y[  =  yi 
P  k 

y2  =  9-—y2-—yi 
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